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Abstract —The concept of All-Back-Schottky-Contact (ABSC) 
thin-film photovoltaic devices is introduced and evaluated using 
2D numerical simulation. The built-in field is generated by reach- 
through Schottky junctions with two metals of different work 
functions in an alternating, side-by-side pattern along the back 
of the device. This approach minimizes the number of interfaces 
through which charge transport is required while eliminating 
the need for transparent conductive front contacts, pn junctions, 
and scribing. Performance metrics are calculated for a range 
of design parameters and it is shown that passivation of the 
front surface and back contacts is critical to achieve efficiencies 
near 20% using only one layer of material deposited on the 
back contacts. This design paradigm may afford a wide range of 
passivation materials since both electrical conduction and optical 
transparency would not be required of any one passivation layer. 

Index Terms —thin-film photovoltaics, numerical simulation, 
semiconductor device modeling 

1. Introduction 

Thin film photovoltaic (TFPV) technologies have the dis¬ 
tinct advantages of high-throughput deposition methods and 
material flexibility. Although TFPV laboratory devices based 
on various materials have recently achieved greater than 20% 
efficiency HI, further advancements are required to tighten 
the gap between cell and module performance, reduce cost, 
and improve long-term reliability. Conventional TFPV cells 
are comprised of several stacked layers of semiconductors and 
dielectrics sandwiched between electrodes, with a transparent 
conducting oxide (TCO) acting as the electrode on the illu¬ 
minated side. Other layers are required to form pn-junctions 
and passivation regions which improve performance. However, 
with each layer comes interfaces which always carry with them 
high densities of recombination centers, energy barriers, and 
seeds of degradation 13 . Performance limiting effects are also 
spawned by the unique physics of inherently nonuniform films 
that are on the order of 1 jim in thickness O. 

As an alternative to the conventional stacked approach, this 
work uses numerical simulation to explore the concept of an 
All-Back-Schottky-Contact (ABSC) TFPV device architecture 
wherein two dissimilar metals of different work functions are 
adhered to the non-illuminated side of the device. Reach- 
through Schottky junctions and the difference in work func¬ 
tions of the two metals establish the built-in field. The basic 
philosophy of this design paradigm is to minimize the number 
of interfaces through which charge transport is required. A 
schematic layout is shown in Fig. [Ha). 

Interdigitated back contact (IBC) crystalline silicon (c- 
Si) solar cells were first conceived im in the 1970’s and 
high efficiency modules are now commercially available 0 , 
0- Several groups continue to study this technology which 
typically employs n-type c-Si wafer material of 100-300 /im 


in thickness with alternating p-type and n-type zones along the 
back to affect charge separation and collection iTTll- lfT^ . Thin- 
film Si devices (2-10 /im thick) have also been developed with 
efficiencies of greater than 10% |[T3ll - lfT5ll . All of these back 
contacted approaches employ pn-j unctions of some type. 

Recently, lateral organic heterojunction cells were devel¬ 
oped with IBCs of alternating metals ifT^ . It was shown that 
such a contacting scheme can in fact produce a built-in field 
that effectively separates charge carriers if the metals are of 
sufficiently different work functions. In that work, the back 
contact materials were Al, Au, and M0O3-coated Au with 
measured work functions of 4.1, 5.1, and 5.6 eV respectively. 
That work also demonstrated the feasibility of using UV- 
lithography to create metal lines down to 5 /im wide separated 
by a channel as narrow as 200 nm. Here we investigate 
how such a contacting scheme may be used to fabricate 
high-efficiency TFPV cells with absorber layers comprised of 
materials similar to poly crystalline CdTe or CIGS. 

The important distinction of ABSC cells from conventional 
TFPV is that the reach-through Schottky junctions can gen¬ 
erate an electric field of E > 10^ V/cm throughout the 
absorber [see Fig. [He)], thereby enabling effective electron- 
hole separation without the need for complex doping steps. 
Some of the main advantages of the ABSC design include: 
(1) lack of TCO, pn-junctions, and window layers, thus re¬ 
ducing parasitic photon absorption and interface issues; (2) no 
scribing requirements for module fabrication thereby reducing 
damage and shadow losses; and (3) opportunities for design 
optimization via optical control and back contact shape/size. 
Several of those points could also help to reduce losses in 
going from cell to module. On the other hand, simplification 
of the semiconductor components leads to added complexity 
in the metal contacting scheme, which is among several 
possible challenges, such as, fabricating a high density of 
Schottky contacts, shunting between contacts, and material 
compatibility between metals of different work functions and 
the semiconductor. 

H. Physical Characteristics 

The basic physics of ABSC-TFPV is similar to that of 
metal-semiconductor-metal (MSM) structures, but with unique 
features due to the device architecture. For example, the elec¬ 
tric potential varies in 2D, such that the standard ID energy 
band diagram is insufficient and the 2D bands follow the 
contours in Fig. Wih), where it can be seen that both transverse 
and lateral electric fields are generated, with the magnitudes 
indicated in Fig. Wic). The calculated conduction band (CB) 
surface is illustrated in Fig. [2] where the same parameter 
values were used as in Fig. [T] The steep CB gradient between 
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Fig. 1. ABSC thin-film device: (a) schematic with semiconductor absorber of thickness h atop alternating layers of metal 1 (Ml), insulator (INS), and 
metal 2 (M2); (b) electric potential distribution in volts relative to vacuum with metal work functions of 0 m i = 4.2 V, 0 m 2 = 5.2 V (other parameters are 
provided in Table |TJ with mMi/2 = ' u ; m 2/2 = Wins = h = 1 /^m. (c) electric field magnitude distribution in V/cm for the parameters of (b). 


the metals (lateral direction) and the relatively small gradient 
between the back and front surfaces (transvers direction) are 
evident in Fig. [21 



Fig. 2. 2D conduction band in eV relative to the Fermi level (E = 0 plane) 
at equilibrium. The back of the device is at = 0. 

Stacked devices are always modeled with the directions of 
light propagation and carrier currents parallel to each other. 
That is not the case in ABSC devices, however, and a 2D 
analysis is required to understand the unique interplay between 
drift and diffusion currents relative to the direction of light 
propagation. There are several characterstic lengths to consider 
(in what follows, we assume a p-type semiconductor and refer 
to the variables shown in Fig. [T]). First, the depletion width, 
Wd of the Schottky contact (which we assume to be ideal) 
at M2 should be within the absorber thickness; h <Wd- The 
depletion width is given by, ini 

where q is the elementary charge, 'ipbi is the built-in poten¬ 


tial, eo is the permittivity of free space, is the relative 
permittivity of the semiconductor, and Na is the acceptor 
concentration. In addition, the semiconductor must be thick 
enough to absorb the majority of incident photons; h > a~^ 
where a is the absorption coefficient of the semiconductor. 
Thus the conditions, 

a~^ < h < Wd, ( 2 ) 

provide the optimal range for the thickness. The depletion 
width also plays a role in determining the lateral field between 
the metals which are separated by the dielectric material of 
width Wins • If is the difference in metal work functions, 
a field on the order of E ^ /S.(j)M/"^ins will be generated 
between the metals when, 

Wins < Wd, (3) 

such that the field is not screened by space charges. 

Another important condition is related to the diffusion 
length of minority carriers, Ln = \/Dr, because they must 
migrate from the quasi-neutral zone above Ml to the high- 
field/depleted region near the insulator and, finally, to M2. 

Therefore, _ 

WMI < 2^iirkT/q, (4) 

where the Einstein relation was used for the diffusion coeffi¬ 
cient, D = jj^kT/q with thermal voltage kT/q, minority carrier 
mobility /i, and r the minority carrier lifetime. The factor of 2 
in Eq. 0 accounts for the requirement of carrier migration for 
half of the contact width. The conditions on the width of metal 
M2 appear to be less stringent since it will depend mostly on 
the majority carrier diffusion length Lp\ hence, 

wm2 ^ 21/ p . (5) 

The electric field distribution above M2 will also affect mi¬ 
nority carrier collection and, therefore, determine how wm 2 
will affect performance. Front surface conditions (charge ac¬ 
cumulation, front surface fields, etc.) will strongly affect the 
field above M2 (discussed further below). 
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III. Numerical Simulation Results 

In what follows, power conversion efficiency ( 77 ), open 
circuit voltage (V^oc), short-circuit current (J^c), and fill factor 
(FF) were determined using numerical simulation for a range 
of parameters including contact widths (wmi and wm 2 ), 
absorber layer thickness (h), surface recombination velocity at 
the back contacts (S'^c) and at the front surface (Sfs). Fixed 
values of Wins = 1 and the difference between metal 

work functions of A0 m = 1 eV were used along with the 
parameter set given in Table [H rather generic for p-type TFPV 
such as CdTe, Ea ciGs, ina, czts, hu, and perovskites. 
II 20 I Shockley-Read-Hall (SRH) recombination was assumed 
in the bulk of the absorber with donor-type defects located at 
mid-gap. Also, ideal Schottky contacts were assumed for this 
initial study. 

Performance metrics were calculated using COMSOL Mul¬ 
tiphysics with a 2D geometry. Simulations were performed 
both in the dark and with exposure to AM1.5G spectrum at 
100 mW/cm^ intensity. The photo-generation rates for this 
study simply followed the Beer Law of exponential photon 
absorption with a direct optical absorption co-efficient ED 
[a{X) = Ay/hu — Eg, with A = 10^ cm“^ photon 

energy hu, and band gap Eg], neglecting reflections at the 
front and back layers. More rigorous optical calculations, such 
as light trapping techniques could be included using the same 
simulation platform. 


TABLE I 

Generic Parameter Values. 


Parameter 

Unit 

Value 

Description 

Eg 

eV 

1.5 

band gap 

X 

eV 

4.0 

electron affinity 

Cr 

[0] 

10 

relative permittivity 

Nc 

cm“^ 

IOI8 

effective CB DOS 

Ny 

cm“^ 

1019 

effective VB DOS 

f-^n 

cm^ s“^ 

100 

electron mobility 


cm^ s“^ 

10 

hole mobility 

Na 

cm“^ 

10i4 

acceptor doping 

Nt 

cm“^ 

1014 

defect concentration 

Et 

eV 

0.75 

defect energy above VB 

O'n 

cm^ 

10-12 

electron cross-section 

<Jp 

cm2 

10-15 

hole cross-section 

(t>Ml 

eV 

4.2 

metal Ml work function 

4>M2 

eV 

5.2 

metal M2 work function 


Note: CB, VB, and DOS stand for conduction band, valence band, and 
density of states, respectively. 


We first consider optimistic surface recombination by set¬ 
ting Sfs = 0 and S^c = 100 cm/s. Bulk SRH recom¬ 
bination is significant with a minority carrier lifetime of 
r = {(TnNtVth)~^ = 1 ns, using the parameters in Table 
U and thermal velocity Vth = 10^ cm/s. As shown in Fig. 
O the efficiency calculations within the parameter space of 
0.5 < h < 4 /im and 1 < wmi = "^ 1^2 < 10 /im 
were within the range of 13% < rj < 18% with the largest 
value given by an optimum thickness of h ^ 1 //m, which 
satisfies Eq. ©; for our parameter values, Wd ^ ^ A^ni and 
a~^ ^0.1 /im. Also, narrower contact width wmi results in 
better performance according to Eq. given our minority 
carrier diffusion length of 0.5 /im. 

Fig. [ 4 ] shows the resulting increase in efficiency when the 


Efficiency (%) 



Contact Width (}im) 

Lig. 3. Efficiency contours as a function of absorber layer thickness and 
contact widths (assuming wmi = '^M 2 ) and the parameters of Table U with 
front surface and back contact recombination velocities of Sfs = 0 and 
Si)c = too cm/s, respectively. Parameters corresponding to the circled area 
at Wmi = wm 2 = 2 are further evaluated in Pigs. [5] and [6| 

minority carrier lifetime is increased to r = 10 ns (by 
lowering the mid-gap defect concentration to At = 10^^ 
cm“^). The efficiency spans from 9% to 21% over a range 
of contact widths 2 < wmi = wm 2 <22 /im and absorber 
thicknesses 0.25 < h < 1.5 /im. In this case, efficiency 
remains above 18% for up to 20-/im contact width when the 
absorber thickness is optimized at < 1. According to Eq. 0 
the allowable contact width increases relative to the case with 
r = 1 ns shown in Fig. [51 The optimum thickness decreased 
slightly below 1 /im because the lower density of occupied 
defect states affects space charge and slightly lowers the built- 
in field which, in turn, lowers ILd in Eq. ©. 

Efficiency (%) 



Contact Width (urn) 

Pig. 4. Efficiency contours as a function of absorber layer thickness and 
contact widths (assuming wmi = wm 2 ) with the same parameters as Fig. [ 5 ] 
except for an increase of minority carrier lifetime from r = lnstOT = 10 
ns. 

It should be noted that Eq. [D may be an oversimplification 
of the depletion width because Wd will also depend on: (i) 
the compensation of Na by charged defect states; and (ii) 
the insulating boundary condition at the front surface which 
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requires that the normal component of the field satisfies n • 
eE = Qt, where n is the surface normal unit vector and Qt is 
the surface charge density due to traps there. As stated below, 
we have assumed neutral surface traps (Qt = 0 ) and, therefore, 
a fiat electric potential at the front boundary which reduces the 
field throughout the entire absorber thickness above M2. The 
properties of both the bulk and surface defects play important 
roles in determining the field distribution E(x). 

Fig. [5] shows the effects of S^c on the performance metrics 
assuming the dimensions of the circled point in Fig. [3] and 
minority carrier lifetimes of r = 1 ns and r = 10 ns. 
Maximum 77 of 19.5% and 21.5% were achieved at S^c = 10^ 
cm/s, which then decayed to nearly stable values between 
15.5% and 16% when S^c = 10^ cm/s (this scenario uses 
Sfs = 0 cm/s). It is apparent that S^c affected efficiency 
mainly through FF and Voc- The latter monotonically decayed 
from 1.0 V to 0.9 V while FF displayed large, non-monotonic 
variations that correlated with 77 . The non-montonic behavior 
of 77 and FF was due to the effect of S^c on the minority 
carrier distribution and, therefore, on the electric field (as is the 
case in a high-injection scenario, n > Na). Our calculations 
show that increasing S^c decreases the transverse field at the 
metal contacts while increasing the lateral field above the 
insulator, resulting in an optimum Sbc- 



Fig. 5. Performance metrics vs. back contact surface recombination velocity 
Sbc assuming S'/s = 0 at the front surface and minority carrier lifetimes of 
r = 1 ns (parameters from Tabled and r = 10 ns. wmi = wm 2 = 2 /xm 
and h = 1 iim. 

Next we consider front surface effects as governed by 
Sfs = (TsVthNst- We set equal capture cross-sections for 
electrons and holes, cFq = 10 “^^ cm^, Vth = 10 ^ cm/s, 
and allowed the surface trap density to vary in the range 
10 ^ > Nst > 10 ^^ cm“^ so that 10 ^ > Sfs > 10 ^ cm/s. 
As an initial test, we assumed only neutral centers located at 


mid-gap. We also considered both the worse case scenario 
of Sbc = 10^ cm/s and passivated S^c = 10^ cm/s. The 
profound effect of front surface recombination is evident in 
Fig. [bias 77 drops from 15.3% to 5.3% over the range of Sfs 
for Sbc = 10^ cm/s. Improving the back surface passivation 
so that Sbc = 10^ cm/s increased 77 to a maximum of 19.1% 
but resulted in no performance gain at large values of Sfs, 
which dominated efficiency reduction. Losses were greatest 
for FF and Jsc, while Vqc exhibited a relatively small loss. 
Qualitatively, those variations are a result of the increase in 
recombination in the area of greatest photo-absorption (at the 
front) while the Schottky junction field is strongest near the 
back. 



Sf3 (cm/s) 

Fig. 6. Performance metrics vs. front surface recombination velocity Sfs for 
Si)c = 10^ cm/s and Si^c = 10^ cm/s at the back contacts. Device parameters 
are from Table U with r = 1 ns, wmi = '^M 2 = 2 iim, and h = 1 fim. 


IV. Discussion 

From Figs. [5] and [ 6 ] it is clear that the challenges of dealing 
with interface quality still remain, but the ABSC cell allows for 
passivation of the front surface by materials that need not carry 
photo-current (i.e. high resistivity and optically transparent). 
Conventional wisdom suggests that a PV cell should have the 
greatest photo-absorption in the region of strongest electric 
field to optimize electron-hole pair separation. Our approach 
goes against that logic but we have shown that a high perfor¬ 
mance device may still be possible due to the presence of both 
lateral and transverse built-in electric fields. An added benefit 
of the ABSC design is that the Schottky junctions are buried 
and therefore protected against the degrading effects of high 
carrier injection. El 

Experimental efforts to develop a proof of concept is of 
foremost importance but additional theoretical analysis may 
include a broader exploration of the parameter space, optimal 


























5 


band gap grading, light trapping measures, low level light 
concentration, and a front surface held layer. Such measures 
could improve efficiency while making the absorber layer 
thinner and the contacts wider. Also, the simulation work 
presented here did not include the resistance of the metal 
contacts, which must be considered when specihc designs are 
studied. 

V. Conclusions 

In summary, an alternative design paradigm of all-back- 
Schottky-contact TFPV has been introduced. The numeri¬ 
cal analysis results suggest that this approach may lead to 
simple yet high efficiency, durable TFPV devices if surface 
recombination can be effectively passivated. In addition, we 
have demonstrated the use of COMSOL Multiphysics® as 
a numerical simulation platform for 2D modeling of TFPV 
devices. 
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